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predicts a non-monotonic dependence of the reduced excess viscosity, vr, on concentration from 
the behaviour of the static structure factor in polyelectrolyte solutions. The theory predicts that 
the peak in r\R occurs at concentrations slightly lower than the overlap threshold concentration, 
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increased concentrations of added salt. The position of the peak, as a function of concentration 
divided by c* is independent of salt concentration or molecular weight. The predictions can be 
tested experimentally. 
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I. INTRODUCTION 

Polyelectrolyte solutions are widely considered to be 
one of the least understood substances in polymer sci- 
ence [ij . There are several features of these solutions that 
make them rather complex. For one, the long-ranged 
nature of the electrostatic interactions results in long- 
ranged correlations even in dilute solutions. In addition, 
polymer conformations are very sensitive to concentra- 
tion and ionic strength because the electrostatic interac- 
tions compete with short-ranged "hydrophobic" interac- 
tions. This complex static behaviour is accompanied by 
very interesting dynamical behaviour of polyelectrolyte 
solutions. In this work, we present a theoretical study of 
the viscosity of dilute and semi-dilute rod-like polyelec- 
trolyte solutions using a liquid state approach. 

The viscosity of polyelectrolyte solutions displays an 
interesting "anomalous" concentration dependence at 
low polymer concentrations 0, ■ The quantity that is 
normally discussed is the reduced viscosity, r)R, defined as 
j]r = (ij — Tfo)/(rfoCp), where r\ is the viscosity of the solu- 
tion, r/o is the viscosity of the solvent in the absence of the 
polymer, and c p is the (monomer) concentration of poly- 
mers. Experiments show that rjn displays a sharp peak at 
low polyelectrolyte concentrations for a variety of differ- 
ent solutions. This anomalous concentration dependence 
of viscosity of dilute polyelectrolyte solutions has been 
the focus of attention for over 50 years and, although 
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there have been many theories that address the problem, 
it is not considered to be well understood 0, 0, • I n 
fact, as discussed eloquently by Forster and Schmidt 1], 
this problem has witnessed some grave mistakes, causing 
much confusion. 

A contributing factor to the above-mentioned confu- 
sion is the fact that a review of experimental data on 
the viscosity of polyelectrolyte solutions shows large in- 
consistencies 0,1H- It is only recently that the reason 
for this discrepancy between different experimental mea- 
surements has been established. In careful measurements 
of the shear rate dependence of the viscosity, Boris and 
Colby 0] and Krause et al. Q showed that polyelec- 
trolyte solutions were shear thinning at extremely low 
shear rates, and argued that most of the older experi- 
ments did not report the relevant viscosity in the low 
shear limit. 

In the last ten years, several theories studies have been 
put forward. Notable among them are the mode cou- 
pling theory calculation by Borsali et al. [(| and by Ra- 
bin et al.Q, the scaling theory of Dobrynin et al. 0], 
the effective medium theory of Muthukumar Q , and the 
Kirkwood theory of Nishida et al. [jj. The scaling the- 
ory does not predict a peak in the reduced viscosity, in 
salt-free solutions. The theory of Muthukumar argues 
that the peak in tjr arises from a screening of intra- 
molecular hydrodynamic interactions as the concentra- 
tion is increased. The theories of Borsali et al. Rabin 
et al.Q, and Nishida et al. predict that the peak in 
r]n arises due to increased screening from countcrions as 
the concentration is increased, and are similar in spirit to 
the mode-coupling theory for charged colloids 10] which 
is argued to be accurate for the viscosity of spherical 
polyelectrolytes Other theories rely on conforma- 
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tional changes of polymers with changing concentration. 
These theories do not take into account the interesting 
behaviour in the static structure factor in dilute polyelec- 
trolyte solutions. 

The basic idea of the present work is that the features 
in rjji arise from the behaviour of the static structure 
factor of dilute and semi-dilute polyelcctrolyte solutions. 
In dilute solutions, the static structure factor displays 
a prominent peak at low wavevectors [T2I llflj . As the 
concentration is increased the peak broadens and moves 
to higher wavevectors (see Fig|21 and discussion in Sec- 
tion III). This indicates the presence of strong liquid like 
correlations on long length scales at low concentrations; 
correlations that become less important as the concentra- 
tion is increased. This observation naturally leads to the 
question: Could this strong non-monotonic concentra- 
tion dependence of static pair correlations be the main 
physics behind the anomalous behaviour of the viscos- 
ity? To this end, we develop a liquid state theory that 
incorporates the behaviour of the static structure fac- 
tor of poly electrolyte solutions. Such a theory, to our 
knowledge, has not been developed previously. While it 
has been suggested [3, HH that the interesting concen- 
tration dependence of viscosity in dilute polyelectrolyte 
solutions could arise from the intermolecular pair corre- 
lations (as reflected in the peak in the structure factor 
at small wavenumbers) , the relationship between the vis- 
cosity and the structure factor is by no means obvious. It 
is of interest, therefore to develop a quantitative theory 
relating the static structure to the viscosity. 

We consider a system of charged rods and present 
a self-consistent mode-coupling theory for the viscosity. 
We choose to study rods in order to focus on inter- 
molecular effects. There have been many theories that 
explain the concentration dependence of the reduced vis- 
cosity on intra-molecular effects. Studying a system with 
rigid molecules allows us to isolate inter-molecular effects 
since conformational changes and intra-molecular inter- 
actions are absent. Although we are not aware of experi- 
ments for the viscosity of rod-like polyelectrolytes, these 
are certainly possible, for example on solutions of tobacco 
mosaic virus (TMV) particles. 

We present a self-consistent mode-coupling theory for 
the viscosity of un-entangled polyelectrolyte solutions. 
Starting with the polymer centre-of-mass density as the 
slow variable, we develop an expression for the polymer 
contribution to the viscosity that is identical to that of 
Geszti 16] . The intermediate scattering function is ob- 
tained from the self-consistent mode coupling theory as in 
the approach of Gotze and co-workers yjj • The centre-of- 
mass structure factor, required in the theory, is obtained 
approximately from integral equation theory [iH, HH |20) . 

The theory explains the behaviour of the viscosity 
based purely on the behaviour of the static structure fac- 
tor. For short chains, the theory predicts that in salt- 
free solutions rjn decreases monotonically with increas- 
ing concentration. With small amounts of added salt rju 
displays a shallow peak as a function of concentration at 



low concentrations. For longer chains, the theory predicts 
a peak in vr as a function of concentration for all salt 
concentrations, including salt-free solutions. The peak 
occurs at concentrations slightly lower than the overlap 
threshold concentration, c*. With the addition of salt, 
the intensity of the peak diminishes, but the position is 
unchanged. For a given salt concentration, the height 
of the peak increases dramatically as the degree of poly- 
merization is increased, but the position is unchanged. 
These predictions can be tested experimentally. 

The rest of this paper is organized as follows: section II 
outlines the theory, section III presents some results and 
a discussion, and section IV presents some conclusions. 



II. MODE COUPLING THEORY 

The polyelectrolyte solution consists of charged rod- 
like polymers and counterions. Each rod consists of N 
tangent charged hard spheres of diameter a; counterions 
are charged hard spheres of diameter a. We combine a 
recently developed quantitatively accurate theory for the 
liquid structure in charged interacting polyelectrolytes 
with a self-consistent mode coupling theory 
(MCT) to study the dynamics. The starting point for the 
calculation of viscosity, 77, is the Green-Kubo formula [2l| 

* = ^ n TJrv f° d< (^(MV^C-M)) , (n.l) 

fc^O fc B J v j Q 

where fee is Boltzmann's constant, T is the tempera- 
ture, V is the volume, a zx (k, t) is the transverse (or off- 
diagonal) component of the wavevector (k) and time (t) 
dependent stress tensor, and (■■•) denotes an average 
over an equilibrium ensemble. The total transverse stress 
tensor of a polyelectrolyte solution contains contributions 
from solvent, polymer, and small ions. 

For dilute and semi-dilute polyelectrolyte solutions 
several simplifications are possible. First of all, the con- 
tribution of the solvent is simply given by the viscosity in 
the absence of the solute, r/ . It is because the presence 
of low concentrations of polymer and electrolyte are not 
expected to alter the solvent dynamics. Secondly, there 
is a contribution of the rotational Brownian motion[22||. 
Lastly, there is a contribution due to polymer-polymer 
interactions rj p - p which is expected to dominate over the 
contributions from small ions. Therefore we argue that 
T) w 770 + rj r + rjp-p and focus on the calculation of the 
polymer contribution. r/ r is calculated by neglecting the 
effect of the interactions on the rotational degree of free- 
dom and is given bv|22| 



where Cr is the rotational friction coefficient and is evalu- 
ated by a hydrodynamic calculation of an ellipsoid of the 
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aspect ratio TV as|2£ 



M(k,t) is the memory kernel given by 



2iri]<j 3 



N 4 



^^(N + VW—i) 
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(II.3) 



For r)p-p, we employ a mode coupling treatment similar 
to that of Geszti [ljj to derive a microscopic expression. 

The first step in the mode coupling approach is the 
choice of the slow collective variables for the description 
of the dynamics of the required correlation functions. 
Natural choice is the hydrodynamic variables, i.e., the 
three momentum current densities of polyion, J a (k) , for 
the co-ordinates a=x,y, and z, and the polyion number 
density, pp(k), defined as pp (k) = J2l=i e ~ tk r * where r i 



M(k,t) 



p P D f dq 



2*0 



F k 2 (q,k-q)f(|k-q|,t)F(g,i), 

(II.6) 

where pp = (pp(k = 0)) = c p /N is the average number 
density of polyion and V (q, k — q) is the vertex function 
given in terms of the direct correlation function c(q) as 

V(q, k - q) = k ■ qc(«) + k • (k - q)c(|k - q|). (II.7) 



is the position of the center of mass of the i polyion 
and Np is the number of polyions. A calculable mi- 
croscopic relation for the viscosity is obtained by using 
the projection-operator formalism to re-write the well- 
known Green-Kubo time correlation function expression 
in terms of P and Q operators 21] . The standard ap- 
proximation in the mode-mode coupling expansion is to 
consider the subspace of various binary products of the 
basic slow variables. Among such binary products, the 
odd ones with respect to time inversion do not contribute 
to the viscosity, and only the even combinations can be 
retained. The two obvious choices of the binary prod- 
uct are the density-density term and the current-current 
term. The current terms are expected to decay much 
faster than the density term, due to the friction with 
the surrounding solvent molecules. Thus, we neglect this 
contribution. Finally, all four-particle correlations are 
approximated as the product of two-particle correlations. 
With the above approximations and simplifications, the 
final expression for the zero frequency viscosity is written 
in terms of the density correlation function of the polyion 
as 

k BT r jm4 f°° (S'{k)} 2 (F(k,t)\ 2 

(II.4) 

where S(k) is the static (center-of-mass) structure factor 
of the polyions, F(k, t) is the corresponding intermediate 
scattering function, S'(k) is the derivative of S(k) with 
respect to k. 

In order to evaluate the viscosity, we need the interme- 
diate scattering function F(k, t) for the polyions and it 
should also be evaluated using MCT. As discussed above, 
we again assume that dynamics of counterions and sol- 
vent molecules is decoupled from dynamics of polyion. 
Then, the equation for F(k,t) is expressed in a closed 
form as0 



Ea. (|II.5|) is a standard MCT equation familiar in the 
in the supercooled liquids and colloids 0, 123. l2q com- 
munity. This equation is a nonlinear integro-differential 
equation which has to be solved self-consistently. The nu- 
merical procedure to solve eq. (|II.5(I following the method 
proposed in Ref.[2J|. 

The bare diffusion coefficient Dq is obtained from the 
value for a long ellipsoid calculated from hydrodynamic 
calculations using stick boundary conditions. For the 
ellipsoid with the aspect ratio of N, the total diffusion 
coefficient is given bv|23| 



Dq = 



Du +2D, 



1 



fc R T 



■In 



(n + \/N 2 - l) 



(II. 



Note that this theory considers only the translational 
motion of the rods, which is assumed to be isotropic. 
We argue that we can neglect the anisotropy in trans- 
lation and its coupling to rotation in the concentra- 
tion regimes we consider. We estimated the contri- 
bution from the anisotropy and coupling to diffusion 
as follows. In the dilute limit, if the interaction be- 
tween polyions is neglected, it is possible to solve the 
rotation-coupled diffusion equation and evaluate F(k,t) 
exactly [22, l2fil 1271 l28l l29| . If the ratio between the paral- 
lel (£>||) and perpendicular (D±) diffusion coefficients is 
not very large, the change in the relaxation rate of F(k, t) 
at short times arising from a coupling with rotational 
diffusion is also small. For example, for Du/Dx — 2, 
the relaxation rate of F(k,t) at short times is changed 
by less than 10 % due to the coupling with the rota- 
tional diffusion. At longer times, t > tr, where tr is 
the rotational relaxation time, F{k,t) is simply given 
by exp[— Dok 2 t] where Do is the average diffusion co- 
efficient defined by D = (Z)|| + 2DjJ/3. Thus, decou- 
pling rotation from translation and assuming the transla- 
tion diffusion is isotropic are reasonable approximations 
in dilute solutions. These approximations become ques- 
tionable when the concentration or rod length becomes 
large when, due to the entanglement effects, the rota- 
tional time increases steeply and anisotropy of the trans- 
lational diffusion will be enhanced. This regime, however, 
is far beyond the scope of the present paper [3(J. 



dt 



(II.5) 

where Dq is the bare collective diffusion coefficient. 



III. STATIC PROPERTIES 

To proceed further we require a model for the polyions 
and a means of calculating the static structure of the 
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polyion centers of mass. In this work the molecules are 
modeled as a collection of interaction sites arranged lin- 
early in a rod-like configuration. Each particle consists 
of N tangent charged hard spheres (or sites) with hard 
sphere diameter cr, which is used as the unit of length 
in this paper. Each sphere carries a negative fractional 
charge, fe, where e is the charge on an electron. The 
effect of solvent and small ions is included into the po- 
tential of interaction between sites on the polyelectrolyte 
molecules. The resulting effective potential, f3u(r), is 
given by 

a ( \ J 00 for r < cr , . 

pu{r) = < „ , v , ( (in.i) 

v ' 11 exp(— nr)/r lor r > cr, v ' 

where = l/k B T, V = fH B /{l + Ka), l B = (3e 2 /e 
is the Bjerrum length, e is the dielectric constant of 
the solvent, and n is the inverse screening length, k — 
y/4:irlB(f 2 c P + 2c s ) where c s is the number density of 
the (monovalent) salt, and c p is the number density of 
polymer sites. In all the calculations presented in this 
work, Ib = 0.758cr, and /=1. If a w 4 A, then an added 
salt concentration of 1 mM corresponds to a reduced salt 
concentration of c s cr 3 s» 4 x 10~ 5 . 

The center of mass static structure factor is calculated 
using integral equations. The single chain structure fac- 
tor, u>(k), is known exactly for this model. The site- 
site static structure factor, S ss (k), is obtained from the 
polymer reference interaction site model (PRISM) the- 
ory |3l|, as described elsewhere ^IJ- ^ nas previously 
been established [33, by direct comparison of theoreti- 
cal predictions for S ss (k) to computer simulations, that 
PRISM is accurate for S ss (k). The center-of-mass struc- 
ture factor is the approximated as S(k) ss S ss (k)/u>(k). 
To check the validity of this approximation, we perform 
Monte Carlo simulations of rods interacting via screened 
Coulomb interactions, and calculate S ss (k), and S(k). 
The simulation algorithm is identical to that described 
elsewhere |32l | except that we do not perform the Ewald 
sum. Figure Q compares simulations results for S(k) 
(filled circles) and S ss (k)/uj(k) (dotted lines) for Ib = cr, 
c p (T 3 = 10 -3 , and 7V=20, and shows that the approxima- 
tion for S(k) is quite accurate. Also shown in the figure 
is the PRISM prediction for S{k). The PRISM S(k) cor- 
rectly reproduces the liquid like structure manifested in 
the peak in S(k). In fact, the theory is in quantitative 
agreement with the simulation results. 

Figure H depicts S(k) from PRISM for iV=150 and for 
several polyion concentrations. We also show the results 
derived from a simpler Debye-Huckel (DH) approxima- 
tion. The DH results is derived by taking the a — > 
limit, and approximating the site-site direct correlation 
function, c ss (r), by c ss (r) = —(3u(r) for all r. The re- 
sulting intermolecular structure factor, denoted S-on{k), 
is given by 

s m (k) J b (in.2) 

1 + WT^ CpU3{k) 
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FIG. 1: Comparison of PRISM predictions (solid lines) for 
the centre-of-mass structure factor to Monte Carlo simulation 
results for c p <t 3 = 10~ 3 , Ib = cr, and iV=20. Dashed lines are 
simulation results for S S s(k)/u)(k). 
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FIG. 2: Static structure factor, S(k), predicted using the 
PRISM theory (solid lines) and the DH approximation (dot- 
ted lines) for N — 150 and c s — ImM, and for various polyion 
concentrations. From left to right, the polyion concentrations 
are c p a [i = 10" 6 , 2 x 10" 5 , 10" 4 , and 5 x 10" 4 . 



Note that, contrary to the PRISM theory, the DH ap- 
proximation is only qualitatively in good agreement at 
very low concentrations, where no strong structure is ob- 
served. 



IV. RESULTS 

For salt-free solutions, the theory predicts that the re- 
duced viscosity tjb. is a monotonically decreasing func- 
tion of polymer concentration, for short chains. As the 
chain length is increased, a peak in tjr is predicted, at 
concentrations slightly below the overlap threshold con- 
centration. Figures|3] (a) and (b) depict r]n as a function 
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FIG. 3: Dependence of the reduced viscosity, r]n, on polymer 
concentration for various salt concentrations, c s — mM (+), 
1 mM (x), 2 mM (*), and 5 mM (□), and for (a) N=20 and 
(b) iV=150. Note that the abscissa is concentration divided 
by the overlap threshold concentration, c* , and both axes are 
logarithmic. 



of polymer concentration for various concentrations of 
added salt, and for degrees of polymerization of 7V=20 
and 150, respectively. (The added salt concentrations of 
1 mM, 2 mM, and 5 mM correspond to reduced salt con- 
centrations of c s a 3 ss 4 xl0~ 5 , 8 xlCT 5 , and 2 xlCT 4 , 
respectively.) In the figures, the abscissa is the polymer 
concentration divided by the overlap threshold concen- 
tration c* which, for this model, is given by c*er 3 = 1/N 2 . 
In all cases, we found that the major contribution comes 
from the polymer-polymer interaction given by the mode- 
coupling expression in eq. I|II.4(I . The contribution of the 
rotational Brownian motion of the individual rod, rj r , is 
independent of the polyion density and does not affect 
the qualitative behavior except for the low concentration 
regime where the mode-coupling contribution becomes 
very small. 

For N=20 (Figure 01 (a)) and c s = (salt-free), t)r is a 
monotonically decreasing function of c p /c*. As the salt 
concentration is increased, the value of rju decreases at all 



polymer concentrations. For iV=20 this results in a shal- 
low peak in r]R at low polymer concentrations and 1 mM 
salt. For high values of added salt rjn is a monotonically 
increasing function of polymer concentration. These re- 
sults are typical of cases when the static structure factor 
does not display a very strong peak at low wave- vectors. 
The influence of the long-range liquid like order on the 
dynamic properties is therefore very weak. The predic- 
tions for short chains are qualitatively similar to other 
theories 0, Q that ignore the effect of static structure on 
the dynamic properties. 

As the degree of polymerization, TV, is increased, the 
theory predicts a prominent peak in rjn that occurs at 
a concentration just below the overlap threshold con- 
centration. The amplitude of this peak increases with 
increasing degree of polymerization and decreases with 
increasing salt concentration. This can be seen in figure 
|3](b) which depicts rjn as a function of c p /c* for iV=150. 
In salt-free solutions, the peak in r]R is very prominent. 
The addition of salt dramatically reduces the height of 
the peak, although a peak is clearly present even for high 
(5 mM) salt concentrations. (Note that both axes are 
plotted a logarithmic scale in order to fit all the curves 
on the same figure.) 

The reduced viscosity is a strong function of chain 
length, in a manner that depends on the salt concen- 
tration. Figures 0] (a) and (b) depict t]r as a function of 
polymer concentration for various values of N and c s = 1 
mM and 2 mM, respectively. In both salt-free and added 
salt solutions, the peak in rju grows with increasing de- 
gree of polymerization, but the position of the peak is 
insensitive to the value of N. In salt-free solutions tjr 
is a strong function of degree of polymerization at low 
values of c p /c* but in added salt solutions rjn is almost 
independent of degree of polymerization for low values of 
Cp/c* . Figure [5] depicts the value of r]R at the maximum 
as a function of degree of polymerization for various salt 
concentrations. The peak values is fitted well by a power 
law except for the salt-free case. The molecular weight 
dependence is very strong, much stronger than what is 
obtained for entangled (neutral) polymer melts. The ex- 
ponent decreases dramatically as the salt concentration 
is increased. 

The physical interpretation of these results is that the 
peak in rju arises from intermolecular correlations be- 
tween the rods. The main physical feature that is input 
into the theory is an accurate estimate of the static struc- 
ture factor of the polyelectrolyte solutions. The viscos- 
ity is then calculated using a fully self-consistent mode- 
coupling theory. Any scaling analysis of the dependence 
of rjR on N and c p must take into account the complex 
dependence of static correlations on the dynamics. 



V. DISCUSSION 

The main ingredients of the theory of this work are 
(i) the use of a fully self-consistent mode-coupling ap- 
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FIG. 4: Dependence of the reduced viscosity, t]r, on polymer 
concentration for various degrees of polymerization, N = 20 
(+), 50 (x), 100 (*), and 150 (□), and for (a) c s =\ mM and 
(b) c 3 —2 mM. Note that the abscissa is concentration divided 
by the overlap threshold concentration, c* , and both axes are 
logarithmic. 



proximation (SCMCT), and (ii) accurate estimates of the 
structure of the solution. It is of interest to determine 
how the actual predictions depend on these two compo- 
nents. We compare the predictions of this work with 
those of related, and simpler, theoretical schemes, as de- 
scribed below. 

(i) Comparison of SCMCT with lowest order MCT 
(LMCT). 

The MCT requires an expression for the intermediate 
scattering function, F(k,t), which we obtain from the 
self-consistent mode coupling equation, ea. (|II.5|l . A sim- 
pler approximation amounts to neglecting the memory 
kernel M(k,t) in ea. l|IL5|) . We refer to this approxima- 
tion, where F(k, t) is simply given by 



F(k,t) = S(k)exp 



D k 2 t 



(V.l) 



as the lowest order MCT (LMCT). Such an approxima- 
tion has been previously investigated by others |(| but 
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FIG. 5: iV-dependence of the peak value of r/u for various salt 
concentrations. Dotted lines are power law fits N" , with v= 
5.8, 3, 2.2, and 1.5, for c s = mM, 1 mM, 2 mM, and 5 mM, 
respectively. 



with the Debye-Huckel (DH) approximation for the static 
structure. The theory of Nishida et al. |2( is closely re- 
lated to the LMCT with an approximate (concentration 
independent) structure factor obtained from numerical 
calculations at zero density. In the LMCT, ca. (|II.4Jl is 
readily integrated over time to give 



kv,T 



?7j?„lmct 



dfcfc' 



[s'(k)y 



(V.2) 



Figure [fj] compares predictions from the SCMCT (solid 
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FIG. 6: Comparison of predictions for the concentration de- 
pendence of rj R from a fully self-consistent MCT (SCMCT) 

(— ) and lowest order MCT (LMCT) ( ) for AT=150. 

PRISM results for S(k) are used in all calculations. In each 
case, the curves correspond to (from top to bottom) salt con- 
centrations of c a = mM (+), 1 mM (x), 2 mM (*), and 5 
mM (□). 

lines) and LMCT (dashed lines) for the concentration de- 
pendence of r/ji for iV=150. PRISM predictions for S(k) 
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are used in both cases. The qualitative behaviour pre- 
dicted by both theories is similar, i.e., there is a peak in 
7]r at concentrations of the order of the overlap threshold 
concentration, and the amplitude of the peak diminishes 
as the concentration of added salt increases. The position 
of this peak is insensitive to concentration of added salt 
and degree of polymerization. The quantitative differ- 
ences, however, are enormous, with SCMCT predicting 
a value of rjji at the peak that is two orders of magni- 
tude higher than LMCT. This emphasizes the fact that 
memory effects arc of considerable importance in these 
systems. This difference between LMCT and SCMCT 
grows with increasing chain length. For example the the- 
ories are almost indistinguishable for iV=20. 
(ii) Influence of static structure. 

In order to see how the accurate estimates of the struc- 
ture affects the results, we compare MCT predictions for 
the viscosity using PRISM results for the structure factor 
to those using the DH approximation for S(k) given by 
eq. (|III.2|I . A good test of the importance of liquid struc- 
ture would be to compare SCMCT with PRISM S{k) 
to SCMCT with the DH S(k). We find, however, that 
SCMCT with the DH S{k) predicts a so-called ergodic- 
nonergodic transition, i.e., F(k 7 t) fails to relax to zero, 
for low polyion concentrations. Such a transition, which 
is also predicted by the MCT for super-cooled liquids, 
leads to a divergence in the viscosity at finite concentra- 
tions! This prediction is clearly incorrect, since no such 
divergence is seen in experiment, or expected on physi- 
cal grounds. This emphasizes, however, that SCMCT is 
very sensitive to the structure factor used as input, as 
one would expect. We attribute the fictitious transition 
to an overestimation of the memory kernel of eq. (|II.6|1 at 
large wavevectors which allows anomalous positive feed- 
back into the relaxation of F(k,t). This over-estimation 
of the memory kernel arises from the broadened and fea- 
tureless S(k) in the DH approximation. In reality, the 
hard-core interaction comes into play at large fc, thus re- 
sulting in a flattening of S(k). The hard-core interaction 
is, of course, neglected in the DH approximation. 

We therefore investigate the influence of structure us- 
ing the LMCT. A combination of the DH and LMCT 
approximations allows us to derive simple scaling results 
for the viscosity. Since the form factor Lu(k) depends on 
k only weakly, we can set oj(k) « oj(0) — N in ea. l|III.2| ) 
to get 



S(k) 



1 



AttIbNc^ 



k 2 



(V.3) 



1 + Aa(x) ' 



where x = k/n, A = AttI b Nc p <j 3 / n 2 = Nc p a 3 /(c p + 2c s ), 
and a(x) = 1/(1 + x 2 ). With these simplifications, t]r 
from the LMCT is given by 



k B TA 2 K 



120ir 2 r] O c p a 3 D 



dxx 



(l + Aa) 2 



(V.4) 



where we have neglected the rotation contribution r] r , 
which does not affect the argument. 

If we ignore the concentration dependence of the inte- 
gral, which is expected to be weak, then 



NWh 

Dn 



-3/2 



(V.5) 



For salt-free solutions, this simple scaling approach pre- 
dicts r]R ~ which is the same as the Fuoss Law or 
the scaling theory of Rubinstein and co-workers 0. In 
solutions with added salt, this scaling approach predicts 
a peak in rj R as a function of c p . The peak occurs for 
c p = 4c s , independent of any of the other parameters. 




FIG. 7: Comparison of predictions for the concentration de- 
pendence of r)n from lowest order MCT using PRISM input 
for S(k) ( — ) and lowest order MCT using Debye-Hiickel input 

for S(k) ( ) for iV=150. In each case, the curves correspond 

to (from top to bottom) salt concentrations of c s =0 mM (+) , 
1 mM (x), 2 mM (*), and 5 mM (□). 



Figure[3compares predictions for r\n using the LMCT, 
and DH and PRISM approximations for S(k). The pre- 
dictions of LMCT with DH S(k) is drastically different 
from that with PRISM S(k). For one, no peak in t]h is 
observed in the DH theory in salt-free solutions. A peak 
appears with added salt, but this peak is very broad com- 
pared to the predictions of the theory with PRISM input. 
The position of the peak moves to smaller concentrations 
as the concentration of added salt is decreased. The pre- 
dictions of LMCT with DH S(q) is qualitatively similar to 
the scaling results, eq.JV3J). The physical reason for the 
peak is completely different from SCMCT with PRISM 
and comes from a competition between the increase in 
rju with increasing concentration and the decrease in rjR 
from decreasing the inverse screening length. The theo- 
ries of Cohen et al. Borsali et al. 0, and Nishida et 
al. are similar in spirit to the lowest order MCT with 
DH input. 
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VI. CONCLUSIONS 

We present a mode-coupling theory for the viscosity of 
dilute polyelectrolyte solutions. Using the static struc- 
ture factor from polyelectrolyte integral equation theory 
as input, we calculate the viscosity using a self-consistent 
mode-coupling approach. The theory predicts a peak in 
the reduced viscosity for long enough polymer chains that 
occurs at concentrations smaller than the overlap thresh- 
old concentration. When the reduced viscosity is plotted 
against the polymer concentration divided by the over- 
lap threshold concentration, the peak position is indepen- 
dent of degree of polymerization or added salt concentra- 
tion, although the height of the peak increases strongly 
with increasing degree of polymerization and decreases 
strongly with increasing the concentration of added salt. 

An important caveat is that this theory does not at- 
tempt to treat entanglement effects. The theory assumes 
that rotational and translational diffusion are de-coupled, 
and that the anisotropy in translational diffusion can be 
ignored. These approximations are not expected to be 
important in dilute solutions but preclude the application 
of this theory to semi-dilute and concentrated solutions 
of polvelectrolvtes [30j . 

The theoretical predictions are sensitive to the level of 
approximation within the MCT and the structure factors 
used. If a lowest order MCT is employed, the theoretical 
predictions are qualitatively similar to the self consistent 
MCT (SCMCT), but the value of the reduced viscosity is 
different by an amount that increases dramatically with 
chain length. If the Debye-Huckel (DH) approximation 
is used for the structure factor, SCMCT predicts a glass 
transition at low concentrations, which arises from the 
neglect of hard-core interactions. With the DH struc- 
ture factor and lowest order MCT (LMCT), the theory 
predicts peaks in the reduced viscosity with polymer con- 
centration, but the position and amplitude of these peaks 
are different from what is observed in the full theory. 

A rigorous comparison with experiment is not pos- 
sible because detailed experiments on the viscosity of 
rod-like polyelectrolytes in good solvents are not avail- 
able. We therefore compare theory to viscosity data of 
(poly)styrene sulfonate (PSS) dissolved in water. Ad- 
ditional complications in these systems are the possible 
importance of intra-molecular effects due to chain flex- 
ibility and hydrophobic interactions. There are signifi- 
cant differences between theoretical predictions for rods 
and experiments on these flexible chains. In many ex- 
periments the peak in ijji occurs in the semi-dilute 
regime, while the focus of theoretical work has been on 
such a peak in the dilute solution regime. The SCMCT 



does predict a second peak in t]r in the semi-dilute regime 
(see for example Figure 01 as well, but, strictly, this is 
outside the regime where we expect the theory to be accu- 
rate. In experiments on sodium PSS |14|. the peak in tjr 
moves to higher concentrations as the molecular weight 
is increased. In the SCMCT for rods the peak position 
occurs at c p = c* for all N, i.e., the peak positions moves 
drastically to lower concentrations as N is increased. We 
carry out the MCT calculation in the same fashion but 
with the center-of-mass structure factor for flexible chains 
and find that it gives essentially similar results as rods 
and does not explain the A-dependence of t^r. This is 
clearly due to the neglect of all other slow modes that 
originate from the intra-molecular configuration appear- 
ing in the case of flexible polymers. Muthukumar § 
makes a similar prediction. The colloid-like theories, in- 
cluding those considered by Cohen et al. Q and Antoni- 
etti et al. [l4|, predict that the peak is independent of 
iV", and in fact occurs for a concentration c p — 4c s , which 
is not in agreement with experiment on PSS. Therefore 
none of the theories can claim to explain experimentally 
observed results, and the occurrence of the peak in T)r 
must still be considered "anomalous" . 

In conclusion, we present a fully self-consistent mode- 
coupling theory for the viscosity of rod-like polyelec- 
trolyte solutions. Significant differences are present be- 
tween the predictions of the theory and experiments for 
flexible chains. Measurements on theoretically well un- 
derstood systems such as solutions of tobacco mosaic 
virus (TMV) will therefore be very useful as a test of 
the theory. Note that since the only input into the the- 
ory is the static structure factor, which has been mea- 
sured in TMV solutions, this theory can be used to make 
parameter-free predictions for the absolute viscosity in 
those systems. We hope this work will stimulate further 
experimental and computational work on this problem. 
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